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Abstract 

Scanning tunneling microscopy (STM) has made tremendous progress in the study and understanding of both classical 
and unconventional superconductors. This has motivated a rapidly growing effort to apply the same techniques to the 
iron-based high-Tc superconductors since their discovery in 2008. Four years have brought exciting advances in 
imaging and spectroscopic investigation of this new class of materials. In this review, we focus on several recent STM 
contributions to the identification of the gap symmetry and the pairing glue. We highlight the unique capabilities and 
challenges still ahead for STM studies of iron-based superconductors. 
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1. Introduction 

Magnetism has long been thought to be antagonistic 
to superconductivity. Therefore, the 2008 discovery 
of superconductivity in iron-containing LaOi_^FvFeAs 
Q, with transition temperature Tc rapidly climbing 
to 55 K upon replacement of La by magnetic rare 
earth elements f2], was unexpected and provoked 
worldwide excitement. As in the case of cuprates, the 
superconductivity emerges from the antiferromagnetic 
parent compounds, suggesting a link between spin 
fluctuations and electron pairing in both materials. 
This finding opens a new avenue to address the 
two ultimate goals in the field of superconductivity: 
to seek the microscopic origin of the high T,. and 
then design new materials with higher r^. Despite 
extensive experimental and theoretical explorations, 
iron-based superconductors (Fe-SCs) still face fierce 
debates on a number of issues including gap symmetry, 
a prerequisite for understanding the secret of high-r^ 
superconductivity. An excellent review by Johnston 
has thoroughly discussed the puzzle in Fe-SCs through 
2010 p]. Another comprehensive review by Stewart 
covers the rapidly-moving field through 201 1 [4|. 

Scanning tunneling microscopy (STM) provides 
unique capabilities to image the atomic and electronic 
structure of a surface with a sub-unit-cell spatial 
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resolution. STM has been applied with great 
success to the study of conventional and cuprate 
superconductors |5 1. The local density of states directly 
measured by STM spectroscopy provides indispensable 
information about the superconducting gap structure, 
its spatial inhomogeneity and behaviors near impurity 
as well as vortex core states. Furthermore, via 
quasiparticle interference (QPI) imaging, STM can 
provide momentum-resolved information about pairing 
symmetry, collective excitations, and competing phases. 
All these accomplishments motivate the growing use of 
STM to Fe-SCs in the quest to understand the intricate 
electron pairing in these materials. 

Following an early review of STM studies of Fe-SCs 
by Yin et al f6l, STM has made considerable new 
progress and greatly contributed to the study of some 
of the most unusual and remarkable properties of 
these materials. A more thorough review by Hoffman 
recently discussed the cleaved surface configurations, 
superconducting and other spectral gaps, and vortex 
states fT\. Our new review here will begin with a 
brief introduction to the STM technique in section 2, 
then concentrate on the experimental highlights in the 
pairing symmetry of Fe-SCs obtained over the past year 
by STM, such as tunneling spectroscopy in section 3, 
quasiparticle interference in section 4, and the vortex 
state in section 5. In section 6, we conclude by 
mentioning a few very recent hints of higher Tc in 
iron-based materials, as well as suggestions for future 
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STM experiments that should shed additional Ught on 
these materials. 

2. Scanning tunneling microscopy and spectroscopy 

STM is based on the quantum tunneling of electrons 
between two electrodes separated by a thin potential 
barrier. A sharp metallic tip, which acts as a local 
probe, is brought within a short distance (typically 
several angstroms) of an electrically conducting sample 
surface. The tip can be positioned with sub-angstrom 
precision in both the xy plane and the z direction 
using a three-dimensional piezoelectric scanner, as 
schematically illustrated in Fig. 1(a). Applying a bias 
voltage between the metallic tip and conducting sample 
leads to a measurable tunneling current; the polarity 
of the bias voltage determines the direction of the net 
electron flow. For instance, a negative bias voltage 
applied to the sample will allow electrons to tunnel from 
the occupied states of the sample through the vacuum 
barrier into the empty states of the tip [Fig. 1(b)]. Upon 
reversing the bias polarity, the electrons will tunnel in 
the opposite direction, from occupied states of the tip 
into empty states of the sample. The tunneling current / 
is well approximated by 



/ oz e-^"'', K = ^^-^ ^ 0.5 V0(eV) A^^ (1) 

where is a mixture of the work functions of the tip 
and sample and d is the tip-sample separation Q. For 
typical metals, (p ~ 5 eV, so / will increase by about one 
order of magnitude for every angstrom decrease in d. 

In a topographic mode, the surface is mapped based 
on the decay of tunneling current / with increasing 
tip-sample separation d. With the bias voltage V fixed, 
the error signal between the measured current /meas and 
the setpoint current /set is fed back to the z piezo to 
control d as the tip is rastered across the sample surface 
[Fig. 1(a)]. The z trajectory of the tip therefore maps a 
contour of constant integrated electron density of states 
(DOS). This technique is referred as constant-current 
mode. In the case of homogeneous DOS, the contour 
corresponds to the geometric topography of the surface. 
However, if the DOS varies spatially, the resulting 
image contains a mixture of DOS and geometric 
information. By choosing the appropriate tunneling 
parameters, the contribution of the inhomogeneous 
DOS can be significantly reduced, and the geometry 
will dominate the image as desired. 

In addition to revealing the geometry of a sample 
surface, STM can also probe the evolution of sample 




Figure 1: (color online) (a) Schematic of an STM. A voltage V 
is applied between a sharp metallic tip and a conducting sample 
surface, leading to a measurable tunneling current /meas which 
decays exponentially with tip-sample separation d. In standard 
constant-current topographic imaging, the difference (or the error 
signal) between /mcas and the setpoint current /jet is fed back to the 
z. piezo to control the tip height, (b) Tunneling process of electrons 
between the tip and sample across a vacuum barrier of width d and 
height (p. The electron wave functions of both the sample and tip decay 
exponentially into vacuum with a small overlap, allowing electrons to 
tunnel between them. With a negative bias voltage V applied to the 
sample, electrons tunnel from the occupied states of the sample to the 
empty states of the tip. 

DOS with energy, up to several eV from the Fermi 
level in both occupied and unoccupied states. The DOS 
can be accessed by switching off the feedback circuit 
to fix d, sweeping the bias voltage V, and recording 
the tunneling current /(V). The conductance dl/dV 
can be obtained either by numerical differentiation of 
I{V) or by a lock-in amplifier technique. In the latter 
case, a small modulation is added to the bias voltage 
V, and the tunneling current / is demodulated to yield 
dl/dV. Although the interpretation of dl/dV spectra 
can be quite complex, in ideal conditions dl/dV is 
a good measure of the sample DOS. If these dl/dV 
spectra are recorded on a dense array of locations in 
real space, the spatial variation of the sample DOS can 
be extracted. This DOS mapping technique has been 
applied to measure local gap variations and magnetic 
vortices (|5]|6llll. 

3. Tunneling spectroscopy 

Like cuprates, all Fe-SCs are layered compounds. 
Most can be mechanically cleaved to expose an 
atomically flat and clean o/j-surface, suitable for 
characterization by surface-sensitive probes. To date, 
STM spectroscopy techniques have revealed a wide 
distribution of reduced gap values lA/k^Tc across 
materials [7J, varying degrees of gap inhomogeneity 
within materials (inhomogeneous doped AFe2As2 IS 
in [lOl vs. homogeneous LiFeAs and FeSe ifTTl 
[12]), multiple superconducting gaps within a single 
material due to the multi-band electronic structure 
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ifTOll . non-universal pairing gap symmetry lfT2l[T3l . and 
symmetry breaking in both superconducting and parent 
phases lfT2l fT4l \T5i . Here, we extend the discussion 
of pairing and normal state symmetry based on recent 
results from higher energy conductance features, atomic 
impurity imaging, and normal state spectroscopy. 



observed in cuprates. However, here Q. - E^ip - A ~ 
(Q.5-2)kBTc (and even the larger iihump - A ~ (2-3)^^7^) 
are much lower than the universal neutron resonance 
energies ~ AAhgT^ found in Fe-SCs ||25l . It is 
possible that surface contamination played a role in this 
study of uncleaved SmFeAsOo.8Fo.2. 



3.1. Electron-boson coupling 

Collective modes, which couple strongly to electrons 
and may serve as "glue" for Cooper pairing, often 
appear as additional conductance features at energies 
beyond the superconducting gap. This higher energy 
structure was first observed and quantitatively analyzed 
in conventional superconductor Pb-insulator-Pb tunnel 
junctions, where minima in dl^/dV^ were found to 
correspond to phonon mode energies lfT6l (It should 
be noted that whether the collective mode energy 
corresponds exactly to a minimum or maximum in 
dl/dV or df/dV^ depends on a detailed calculation 
which must be carried out specific to each material). 
Subsequent tunneling experiments on cuprates showed 
a similar dip-hump structure in dl/dV at energies above 
the gap (5], and several studies specifically matched 
the spectral features I17l tl8J to the energy of a spin 
resonance mode D.r near (n, n) observed by neutron 
scattering [T9l l20ll . This supports the importance of 
spin fluctuations in the pairing mechanism of high-Tc 
cuprates. 

Soon after the discovery of Fe-SCs, the importance 
of spin fluctuations was suggested in these materials 
as well, with the proposal that the high-r^ 
superconductivity arises through spin flip quasiparticle 
excitations between the F-centered hole pockets and 
M-centered electron pockets (so-called s-t wave pairing) 
ETl I22I . In this context, tunneling spectroscopy is 
again a powerful method to measure the energies of 
collective modes coupling to electrons, and to compare 
these energies to spin resonance modes detected by 
other techniques. 

Fasano et al first observed a dip-hump feature 
by spectroscopic STM in nearly optimally doped 
1111-type SmFeAsOo.8Fo.2 {T,- =45 K) [Fig. 2(a)], and 
argued that this feature was likely a spin resonance 
mode ||23]| . Unlike the superconducting gap A ^ 
7 meV, the outer dip-hump structure is asymmetric 
with respect to the Fermi level {E^). On the positive 
energy side (empty states), the energy of the presumed 
collective mode, Q = isdip - A, ranges from ~2-8 meV, 
and locally anti-correlates with the presumed pairing 
strength A. This phenomenology is similar to the 
local anti-correlation of the presumed phonon [24J and 
spin resonance IfTSl mode energies with gap energies 
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Figure 2: Large-energy-scale tunneling conductance spectra showing 
the dip-hump structure beyond the superconducting gaps in various 
Fe-SCs: (a) SmFeAsOo.gFo.a {T, =45 K) I23j; (b) Bao.6Ko.4Fe2As2 
(T, =38 K) [26]; (c) LiFeAs (T, =17 K) ED; (d) NaFeo.975Coo.o25As 
(Tc =21 K) 123; (e)FeSeo4Teo6 (^f =14.5 K) fl3|; (f) FeSe (T,=93 
K)(n|. 

Shan et al subsequently reported a spectroscopic 
STM study of hole-doped 122-type Bao.6Ko.4Fe2As2 {T^ 
=38 K) 112611, in which features of a collective mode 
were more clearly visible beyond the superconducting 
gap in raw dl/dV spectra below [Fig. 2(b)]. The 
mode energy of Q. » 13-14 meV, calculated as the 
difference between the dP/dV^ minimum and the larger 
superconducting gap A x 8.4 meV, closely agrees with 
the spin resonance excitation Q,. x 14 meV as measured 
by neutron scattering [28 J. Again, Q anti-correlates 
locally with A within this material, which contrasts with 
the global positive correlation between spin resonance 
excitation D.r and T^. across 6 different Fe-SCs [i25l . 
Balatsky et al addressed this same discrepancy in 
cuprates, proposing a local strong coupling model, 
where A(r) not only scales linearly with the boson mode 
energy f2(r), but also correlates exponentially with the 
boson coupling constant g{r) via A(r) ~ Q(r)e"'^*'^''-' 
[29]. In other words, g(r) is more important than Q(r) 
in determining A(r). In addition, it was found that 
the effective coupling constant is inversely proportional 
to the boson energy Q.{r). Therefore, larger Q.(r) 
corresponds to weaker coupling g(r), and consequently 



3 




to smaller A(r). 

To round out the four major families of Fe-SCs, Chi 
et al and Wang et al observed the dip-hump feature 
in 111 -type stoichiometric LiFeAs {T^ -\1 K) Il27l 
and slightly Co-doped NaFeo gvsCoo.oasAs {Tc -21 K) 
single crystals |25|, as illustrated in Figs. 2(c) and 
2(d), respectively. Both the deduced bosonic mode 
energies, Q = 5 meV in LiFeAs and Q = 7.8 meV 
in NaFeo.975Coo.o25As, are close to the spin resonance 
excitation Qr ~ 8 meV QUI . Moreover, such features 
are also visible in dl/dV spectroscopy on 11 -type iron 
chalcogenides FeSeo.4Teo.6 (T,. =14.5 K) [13| and FeSe 
iTc=9.3 K)^ [Figs. 2(e) and 2(f)], although not 
clearly explained in these papers. 
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Figure 3: Collective resonance mode energy SI from STM (filled 
symbols) and spin resonance n, from neutron scattering (open 
symbols) in both cuprates and Fe-SCs vs. bulk T^. The n data 
for Bi2Sr2CaCu208+a comes from dl/dV spectra from 17 dift'erent 
SIS break junctions at various hole dopings, and is computed 
as n = £(dip in rf//dy) - 2A (17) ■ All other Q. data comes 
from STM junctions with normal state tips. Pro 88LaCeo.i2Cu04: 
a = E(peakmdf/dV^) - A (H; SmFeAsOo.gFo.i: H = 
£(dip in dI/dV)-A |23|; Bao.6Ko.4Fe2As2: CI = ^(dip in df/dV^)- 
A HI; LiFeAs: = £(dip in dI/dV)-A (27l; Na(Feo.975Coo.o25)As: 
n = Eidipindf/dV^) - A (25); FeSeo.4Teo.6 and FeSe: £1 = 
ii(peak in dl^/dV^) — A I12II131 . Neutron scattering data comes from 
Refs. (I9]|20]|25]|^|30l. 

Figure 3 summarizes the bosonic mode energy 
Q derived from STM measurements and Q.,- from 
neutron scattering in both cuprates and Fe-SCs, vs. the 
superconducting transition temperature T^.. Although 
changes by a full order of magnitude, it is remarkable 
that most data collapse onto the universal relation Q 
(or D.r)/kBTc = 4.8 ± 0.3. This indicates that the 
dip-hump structure may be of the same origin in both 
material families. In conclusion, the dip-hump structure 
has been resolved by STM in the four major families 
of Fe-SCs and in two families of cuprates. These 
features link to the spin resonance Q.,-. The universal 
relation Q ~ A.MbTc suggests that the interactions of 
electrons with the spin resonance mode are crucial for 



the superconductivity in Fe-SCs. 

3.2. Atomic impurities 

Impurities in superconductors are an active subject 
of study, with the potential to unravel the pairing 
symmetry of the superconducting state. Impurities 
are particularly important in cuprates and Fe-SCs, 
most of which require chemical substitutions to 
enable superconductivity. While in conventional 
s-wave superconductors, only magnetic impurities 
cause pair-breaking and suppress superconductivity, 
in t/-wave superconductors (e.g. cuprates) [31 1 and 
multi-band s± wave superconductors (e.g. possibly 
Fe-SCs) 13211331 [341. both magnetic and non-magnetic 
impurities are expected to induce bound states in the 
superconducting gap. 

In Fe-SCs it has been theoretically suggested that 
impurity studies can address the ongoing debate 
between the s-t and s++ pairing symmetry, since only 
the sign change of the order parameter in s^ pairing can 
give rise to sub-gap bound states around non-magnetic 
impurities ||32|[33l|34l. Therefore, spectroscopic study 
of the energetic structure of impurity-induced bound 
state resonances could constrain the pairing symmetry. 
Furthermore, STM study of the spatial structure of 
impurity-induced states could address proposed phases 
such as orbital ordering which breaks the underlying 
crystal lattice symmetry ll35l l36ll . Exploration of 
single impurities is most promising in nominally 
stoichiometric superconductors with unreconstructed 
surfaces such as FeSe(OOl), KFe2Se2(l 10), and 
LiFeAs(OOl) 01371138]. 

Single atom impurities were first identified on 
stoichiometric FeSe films grown by molecular beam 
epitaxy (MBE) [12|. Both Fe adatoms and Se vacancies 
showed sub-gap resonances in dl/dV at symmetric 
energies, with more pronounced filled state peaks [Figs. 
4(a)-4(c)]. These sub-gap resonances would indeed be 
expected at impurity sites for unconventional s± wave 
pairing in FeSe |[32l l33l l34l . However, the existence 
of resonances here does not prove the s±- pairing, as 
the magnetic nature of the impurities is unknown: if 
magnetic, they would possibly induce bound states even 
in the case of simple i-wave pairing [31 ]. 

Subsequently, Li et al introduced two types 
of mirror-symmetric Fe vacancies [Fig. 4(d)] into 
superconducting KFe2Se2(110) films [37 1. The dl/dV 
spectra on both impurities reveal strongly suppressed 
coherence peaks and a pair of in-gap resonances peaked 
at +1.9 mV [Figs. 4(e) and 4(f)]. Although the energies 
of the resonances at both sites are again symmetric with 
respect to Ep, their strengths are again quite different. 
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Figure 4: Impurity-induced bound states in the superconducting gaps 
of MBE-grown FeSe(OOl) and KFe2Se2(110) films, (a-c) Schematic 
atomic structure of Fe adatom and Se vacancy, and their effects on the 
dl/dV spectra at 0.4 K. Insets show the corresponding topographic 
images (3 nm X 3 nm, Vs = 10 mV, / = 0.1 nA) with white dots 
indicating the topmost Se atoms. Pronounced bound states were 
observed at -1.4 meV and -0.4 meV on Fe adatoms, and at -1.0 meV 
on Se vacancies. The tunnel junction was set at 10 mV and 0.1 nA 
1121 . (d-f) Atomic structure and magnetic-field dependence of dl/dV 
spectra measured on two inequivalent Fe vacancies at 0.4 K. Insets 
show the corresponding topographic images (3 nm X 3 nm, Vs = 30 
mV, I = 0.03 nA) with black dots indicating the K atoms. Black dashes 
depict the evolution of in-gap bound states with the applied magnetic 
field. The tunnel junction was set at 15 mV and 0.1 nA 1371 . 



but with more pronounced empty state peaks in contrast 
to FeSe [12 J. More interestingly, the energies of the 
in-gap bound states on the two vacancies were found 
to shift linearly in opposite directions with applied 
magnetic field [Figs. 4(e) and 4(f)], demonstrating 
that the two Fe vacancies carry spins of different 
orientations. Therefore, the observed in-gap resonances 
do not contradict the theoretical predictions for either Sj- 
or pairing symmetry. 

Very recently, Grothe et al systematically 
characterized five predominant defects in nominally 
stoichiometric LiFeAs single crystals ["381. Despite 
their unknown chemical nature, these defects can be 
named according to their two-dimensional point group 
symmetries, and categorized into two distinct sets. 
One set consists of two different Fe-D2 defects at the 
Fe sites, which preserve the local lattice symmetry, 
and exhibit only a single resonance near the edge 
of the smaller gap in this multi-gap superconductor 
The other set, including Fe-C2, As-Di and Li-Di, 
was found to break the local lattice symmetry, and 
induce additional in-gap bound states, pronounced at 
either positive or negative biases. All five of these 
defects suppress the superconducting coherence peaks, 
and their in-gap resonances show such pronounced 
particle-hole asymmetry in the spectral weight that for 
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Figure 5: Symmetry breaking in Fe-SCs. (a, b) Bound states maps in 
the vicinity of single Fe adatom (1.5 nm X 1.5 nm) and Se vacancy (3 
nm X 3 nm) of FeSe films 1 12|. The white dots indicate the topmost 
Se atoms, (c, d) Geometric (2 nm X 2 nm, = 6 mV, / = 0. 1 
nA) and electronic (6 nm X 8 nm, Vj = 10 mV, / = 0.1 nA) dimers 
induced by the Se substitution for Fe sites in FeSe films |39|. (e) 
Geometric dimers in LiFeAs single crystals (19 nm X 19 nm, Vg = 
20 mV, / = 10 pA) lack accompanying electronic dimers [ll]. (f) 
Average electronic dimer induced by Co substitution at the Fe site in 
Ca(Feo.97Coo.o3)2As2 (5 nm X 5 nm) [40). 

some resonances no corresponding symmetric peak can 
be detected on the opposite side of Ep. Although the 
magnetic nature of the impurities is unknown, one may 
speculate that some of the five may be non-magnetic, 
in which case the consistent existence of sub-gap 
resonances at all five different impurities would match 
theoretical predictions for s±- but not s++ pairing 
symmetry 02l[33ll34ll. 

In addition to the bound state spectra of these 
impurities, which provide important information about 
pairing symmetry, their spatial structure, which often 
breaks the local lattice symmetry, supports the proposed 
orbital ordering phase [35. .36J . In FeSe, taking into 
account the full three-dimensional structure of each 
FeSe sheet, both Fe adatoms and Se vacancies have 
structural C4,> symmetry (aside from the unresolvably 
small 0.4% orthorhombicity lfT2l '). yet both of their 
bound state dl/dV maps show pronounced C21. 
symmetry [Figs. 5(a) and 3(b)]. This likely originates 
from an orbital ordering effect [12]. Furthermore, 
Se substitutions at the Fe site were also identified 
II39I . Although these defects produce stochastically 
oriented short-range "geometric dimer" signatures in 
keeping with the local C2,, symmetry of their lattice 
sites [Fig. 5(c)], they are also found to produce longer 
range a-axis-oriented electronic dimers separated by 
~ 16flFe-Fe [Fig. 5(d)]. Therefore, analogous to Se 
substitutions at the Fe sites in FeSe films [39|, the 
two orthogonal Fe-D2 "geometric dimers" in LiFeAs 
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most probably originate from the two inequivalent 
Fe positions in LiFeAs |38|. Interestingly, the 
Fe-D2 impurities in LiFeAs do not appear to produce 
electronic dimers [Fig. 5(e)] IfTTl . suggesting that 
neither orbital ordering ll35l |^ nor pocket density 
wave ordering, both of which may give rise to such 
dimers BTl . arises in the tetragonal phase of Fe-SCs 
such as LiFeAs. Moreover, although the density of 
Co dopants substituted at the Fe sites in orthorhombic 
Ca(Feo.97Coo.()3)2As2 is too high to clearly distinguish 
spectial feature around individual Co aotms, an 
atomically registered averaging analysis shows that 
each dopant also forms an a-axis-oriented electronic 
dimer [Fig. 5(f)], similar to electronic dimers in FeSe, 
but with peaks separated by 22 A, approximately 

8flFe-Fe IEQI. 

3.3. Normal state 

STM spectroscopy was also used to probe the 
electronic structure of iron-based compounds in their 
parent and normal states. Most such studies have 
revealed a quite broad and asymmetric V-shaped DOS 
suppression around /ip [7], which is reminiscent of 
the controversial pseudogap in cuprates |5|. In 
the case of iron-based materials, however, several 
other phenomena, such as surface reconstructions, or 
structural and magnetic transitions above Tc, may be 
responsible for the DOS suppression around Ep- It 
therefore appears more challenging to understand the 
DOS suppression in iron-based materials. 

Zhou et al recently investigated the doping and 
temperature evolution of the electronic structure of 
Ill-type NaFci-^COfAs single crystals by STM 
spectroscopy [Fig. 6] ll42l . In the parent compound 
NaFeAs, temperature-dependent dl/dV spectra directly 
reveal a gap opening just below the spin density 
wave (SDW) transition temperature Tsdw ~ 40 K. 
The observed gap edge peaks are consistent with the 
gap opening associated with a new electronic ordered 
state (such as superconductivity or SDW), where the 
DOS conservation will push states out toward the gap 
edge. In the case of the un-reconstructed surface of 
cleaved NaFeAs crystals, it is reasonably argued that 
the observed gap correlates closely with the SDW in the 
parent NaFeAs compound. Moreover, the SDW-related 
gap microscopically coexists and competes with the 
superconducting gap in an underdoped NaFei_ jCo^ As 
(x=0.014) sample, which lies at the SDW and 
superconducting phase boundary 1431 . 

With increasing Co doping, the SDW gap is 
suppressed and an inhomogeneous superconducting gap 
occurs in optimally doped NaFei_ jCO j As (x=0.028). In 



contrast to cuprates, no pseudogap-like feature above 
Tc is found in the dl/dV spectra, similar to the doped 
FeSe().4Teo.6(001) [13] and stoichiometric FeSe(OOl), 
LiFeAs(OOl), and KFe2Se2(110) surfaces OH [TTl 
|27]p7|. Unexpectedly, in overdoped NaFei_vCOj:As 
(x=0.061 and 0.075), a spatially inhomogeneous 
pseudogap-like feature re-enters and persists to the 
relatively high temperature of 50 K. Since no gap edge 
peaks can be clearly observed at positive bias voltages 
up to 50 meV, it is not sufficiently convincing to claim 
that the DOS depression around Ep originates from 
yet another ordered state. The nature of the DOS 
depression in overdoped NaFei_jCo^As needs further 
investigation. 
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Figure 6: (a) Temperature dependence of dl/dV spectra in 
NaFei-vCoj As single crystals with various Co concentrations x. The 
blue (red) curve marks SDW (superconducting) transition. The curves 
are vertically offset for clarity |42|. (b, c) V2 X V2 charge ordering 
(Vs = 40 mV, / = 0.02 nA) and dl/dV spectrum in KFe2Se2 (001) films 

Ell 

Using STM, Li et al imaged four different phases in 
MBE-grown KFe2Se2(001) films on SrTiOa substrates. 
They argued that the parent compound corresponds to 
stoichiometric KFe2Se2 with V2 x V2 reconstruction 
with respect to the original Se lattice [Fig. 6(b)] [44|, 
although this surface had previously been suggested 
to be superconducting due to the observed V-shaped 
gap by Cai et al [45J. However, the nonvanishing 
DOS(£'f) at 0.4 K [Fig. 6(c)] suggests that the V-shaped 
gap of Cai et al may have a non-superconducting 
origin Il44l . Moreover, Cai et al explained the V2 x 
V2 reconstruction as the consequence of the block 
antiferromagnetic state of the underlying Fe layer ll45l . 
Combining both studies, the V-shaped gap in the 
reconstructed V2 x V2 surface possibly originates from 
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the SDW transition, similar to NaFeAs Il42ll43]| . 

4. Quasiparticle interference 

In addition to inducing bound states, impurities can 
also scatter quasiparticles and lead to energy-dependent 
QPI patterns, which can be imaged in real space 
with spectroscopic STM |5, 7|. These patterns are 
then Fourier transformed to quantify the periodicity 
of electronic modulations. The combination of the 
real-space QPI pattern and momentum-space analysis 
is often called Fourier-transform scanning tunneling 
spectroscopy (FT-STS). Applications of FT-STS to 
cuprates have revealed dispersive or static electronic 
modulations in vortex cores, superconducting, and 
pseudogap states Q. These findings have shed light on 
the microscopic nature of both the superconducting and 
pseudogap states. FT-STS can even be used to reveal 
the momentum-dependent coherence factors associated 
with quasiparticle scattering and thus distinguish 
different pairing symmetries in unconventional 
superconductors IIT3ll46l . 

In Fe-SCs, FT-STS may present richer QPI patterns 
due to the multiple Fermi surfaces derived from the five 
Fe t/-orbitals |3|. FT-STS on non-superconducting 
Ca(Feo.97Cooo3)2As2 gives evidence for static, 
unidirectional electronic nanostructures aligned 
along the crystalline a axis lfT4l . which are further 
suggested to arise from the combined effects of all 
electronic dimers induced by Co substitutions at Fe 
sites [Fig. 5(f)] flOl . In fact. Song et al more directly 
observed similar electronic dimers induced by Se 
substitutions at Fe sites in superconducting FeSe films 
[Figs. 5(d)] lfT2l l39l . These findings provide evidence 
for the existence of a more complex electronic nematic 
state in the orthorhombic phase of some Fe-SCs, which 
may be responsible for the transport anisotropy in 
these materials ['40', "471. Anisotropy in Fe-SCs without 
Fe-site dopants remains to be understood. 

Recently, two groups performed FT-STS studies of 
LiFeAs, because of its clean and charge neutral cleaved 
surface [48, i49J. Hanke et al first found that the 
Fourier transformed QPI images closely resemble the 
LiFeAs Fermi surface itself [Figs. 7(a) and 7(b)]. The 
QPI images are then explained as a series of interband 
scattering processes between a small F-centered pocket 
(which acts as a van Hove singularity) and other 
pockets 1,48.1 . This contrasts with most other Fe-SCs 
where the scattering between the well-matched hole and 
electron pockets plays a crucial role in both QPI patterns 
and i-t-wave pairing [13|. Furthermore, Hanke et al 
calculated the expected QPI in the superconducting state 



with various elementary order parameters (.?++, d 
and p waves) within the BCS model, and concluded 
that their data supported p wave pairing symmetry in 
LiFeAs, although a more complex order parameter such 
as i -H id wave could not be excluded. It should also be 
noted that the observation of a spin resonance in LiFeAs 
by neutron scattering does not support spin-triplet 
p-wave superconductivity in LiFeAs but instead suggest 
that the mechanism of superconductivity is similar to 
that in the other Fe-SCs BUI . 




Figure 7: (a) Schematic constant energy contours (CECs) of LiFeAs at 
-11.7 meV with the dashes indicating the first Brillouin zone. There 
are two hole-like pockets around F and two electron-like pockets at 
the zone boundary. As marked, q\,,A represent various scattering 
processes among the pockets, while q^ f, represent umklapp processes, 
(b) Fourier transformed image of the measured QPI at the same energy 
as (a). The QPI features match well with the scattering g vectors, 
(c-g) Calculated QPI (q space) in the normal and superconducting 
states with .v+, d, p, and i++ wave pairing symmetries |48|. (h) 
Anisotropic energy gap structure on the three hole-like bands 
hi, ho, and of LiFeAs single crystal at T = 1.2 K. The 0.35 
meV error represents the thermal resolution of FT-STS at 1.2 K. (i) 
Three-dimensional image of the measured A^. (solid dots) on the three 
hole-like bands |49|. Four-fold symmetiic superconducting gaps are 
clearly visible with their minima aligned along FM (FX) for (hi 
and hi) band. 

Subsequently, Allan et al measured intraband QPI in 
LiFeAs to characterize the anisotropic superconducting 
energy gap P9| . They observed three distinct 
holelike F-centered bands hi, /12, and h^,. At specific 
energies, the Fourier transformed images of Bogoliubov 
QPI giq, E) exhibit disconnected slices with vanishing 
scattering intensity in certain directions for each band, 
evidencing anisotropic A^. Analysis shows that the 
gap minima correlate with the locus of maximum 
intensity in g{q,E). Therefore, a plot of intensity 
maxima in a g(q,E) plane along a high symmetry 
direction contains the information on A;;. Figures 7(c) 
and 7(d) summarize the experimental results. The 
superconducting energy gap A^ for the /13 (hi and h2) 
band shows four-fold symmetry with its minima aligned 
along FM (FX) orientation. This differs from FeSe 
where two-fold pairing symmetry was suggested based 
on the anisotropic vortex structure |[T2| . 
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5. Vortex State 

Vortices, quantized tubes of magnetic flux, form as a 
magnetic field penetrates into a type-II superconductor. 
Since the pioneering work by Hess et al \5U\, STM 
has served as a powerful technique to visualize vortex 
arrangements in superconductors @ |6l |71. More 
significantly, STM can uniquely access the internal 
structure of a single vortex core, providing a measure 
of the superconducting coherence length f lH (TO) and 
even fundamental information about electron pairing 

miiini. 

In 1964, CaroU et al theoretically predicted a series of 
quasiparticle bound states separated by ~ A^/Ep within 
the vortex cores iBTI . Since Ep is typically significantly 
larger than A in conventional superconductors ll52ll . 
the excitation spectra are quasi-continuous, and will 
appear as a pronounced peak at Ep, which can split 
into two symmetric peaks and eventually merge into 
the coherence peaks away from vortex center. This 
has been experimentally verified in NbSe2 iBOl |53]| . 
In Fe-SCs, however, due to small carrier density 
(small Ep) and high T, (large A) |3|, A^/Ep is 
relatively large and the excitation spectra are no longer 
continuous. Therefore, the quantum limit regime, where 
the DOS spectrum at the vortex center exhibit strong 
particle-hole asymmetry, may be anticipated Ii52ll54ll . 

STM has revealed bound states within the vortex 
cores of stoichiometric FeSe (TT\ and LiFeAs fTTI, as 
well as hole-doped Bao.6Ko.4Fe2As2 [10|. In FeSe, 
Song et al observed a vortex core state exactly at 
Ep [Fig. 8(a)], which split into two symmetric peaks 
and eventually merged into the coherence peaks along 
the Fe-Fe orthorhombic b axis, analogous to the 
conventional superconductor NbSe2 |53 |. The splitting 
occurs along the Fe-Fe orthorhombic a axis as well, 
but the two symmetric peaks never merge into the 
coherence peaks, providing evidence of the two-fold 
symmetry in electron pairing. In contrast, STM studies 
of LiFeAs and Bao.6Ko.4Fe2As2 revealed significant 
particle-hole asymmetry in the DOS spectrum of the 
bound states, with pronounced peaks below Ep [Figs. 
8(b) and 8(c)]. This suggests that LiFeAs and 
Ba{)(,K()4Fe2As2 may be located in the quantum limit 
regime. 

In addition, STM conductance maps reveal 
material-dependent vortex core shapes. In FeSe, 
the vortex appears elongated along the a-axis [Fig. 
8(d)], leading to the suggestion that pairing is two-fold 
symmetric in FeSe lfT2]| . In contrast, Hanaguri et al 
identified a four-fold star-shaped vortex in LiFeAs, 
with high LDOS tails oriented at 45° to the nearest 



Fe-Fe a and b axes [Fig. 8(e)] ifTTI . A model of 
gap anisotropy fails to account for the star-shaped 
vortex, because the gap minima in the outermost hole 
cylinder governing the LDOS distribution occur along 
the Fe-Fe directions 1 49] \55l and thus should give 
rise to high LDOS tails along a and b axes, in stark 
contrast to the experiments. To address this confusion, 
Wang et al theoretically calculated the shape of vortex 
cores by considering the Fermi surface anisotropy 
ll56l . The results highlighted that the Fermi surface 
anisotropy can dominate the vortex core shape and 
prevent direct access of superconducting gap features 
if the gap itself is not highly anisotropic. Therefore, 
the observed four-fold star vortex core shape could 
be explained solely by the Fermi surface anisotropy 
in LiFeAs. Indeed, the ARPES experiment showed 
that the outermost hole cylinder exhibits a rounded 
square cross section with flat regions at 45° to the Fe-Fe 
directions in LiFeAs ll55l . which should lead to more 
quasiparticles and LDOS tails along these diagonals. 




Figure 8: Vortex bound states and structure, (a-c) dl/dV spectra 
taken at the vortex center in (a) FeSe 1121 . (b) LiFeAs 1111 . (c) 
Bao.6Ko.4Fe2As2 (T0|. (d-f) Vortex structure in (d) FeSe fWi, (e) 
LiFeAs (f) Sro.75Ko.25Fe2As2 j57J. White arrows mark the 

nearest Fe-Fe directions. 

In nonstoichiometric Fe-SCs such as 
BaFe1.8Coo.2As2 181, Bao.6Ko.4Fe2As2 M HD, 
NaFeo.975Coo.o25As 1125] and Sro.75Ko.25Fe2As2 lE), 
single vortices often appear irregularly shaped and 
it is challenging to identify the vortex structure. To 
overcome this issue. Song et al registered the vortex 
centers and then averaged the density of states from 
48 vortices in Sro.75Ko.25Fe2As2 (T^^dl K) (51], 
as depicted in Fig. 8(f). Although the underdoped 
Sro.75Ko.25Fe2As2 lies in the orthorhombic phase Q, 
the averaged vortex appears nearly isotropic, in contrast 
to the two- and four-fold symmetric vortices in FeSe 
and LiFeAs IfTTI [T2l . This indeed supports the claim 
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that the small orthorhombic lattice distortion can't bear 
the responsibility for the anisotropic vortices in FeSe 
IIT2I . Further experiments are needed to understand the 
internal structure of single vortices in doped Fe-SCs 
and their bearing on the pairing symmetry. 

6. Conclusions 

We have shown that in the past year the STM 
has made significant contributions in the quest for 
fundamental understanding of the microscopic electron 
pairing mechanism in Fe-SCs. Low temperature 
tunneling spectroscopy has revealed: (i) a generic 
dip-hump spectral feature beyond the superconducting 
gap in all four main structural families of Fe-SCs lfT2l 
[T3I l23l I25I l26l l27l . which links to the spin resonance 
detected by neutron scattering f25] and supports the 
importance of spin fluctuations in electron pairing in 
Fe-SCs; (ii) sub-gap bound states at single impurities in 
stoichiometric FeSe, LiFeAs and KFezSez Ifnil37ll38l . 
consistent with a s± wave electron pairing picture; 
(iii) rotational symmetry breaking in both parent and 
superconducting states lfT2l fl4l [TSl l40ll . supporting 
an orbital ordering phase in Fe-SCs 1351 |36J ; (iv) 
SDW-related gaps opening near in the Fe-SC parent 
compounds EH |43] EH |45l. In addition to these 
achievements, QPI imaging has provided information 
about superconducting pairing symmetry lfT3l l48l l49l . 
a prerequisite for the final determination of the nature 
of the high-rc superconductivity in Fe-SCs. Finally, 
STM has directly imaged material-dependent vortex 
core states and shapes IfTOlfTTl [1211571 . Such results hint 
at the versatility in electron pairing, and can ultimately 
provide material-specific tests of the pairing symmetry 
in Fe-SCs. 

Following these advances, there are a number of 
remaining questions to be addressed by the STM 
technique. In contrast to cuprates, the multi-band 
nature of Fe-SCs leads to electron pairing which relies 
sensitively on material and doping [3 1. Even within the 
same material, the superconducting order parameter can 
change from nodal, in thick FeSe films on graphene 
lfT2ll . to nodeless, in single unit-cell FeSe films on 
SrTiOs p58l . The grandest challenge remains to clarify 
the pairing symmetry and its origin, as a prerequisite to 
understanding the secret of high-Tc superconductivity 
in Fe-SCs. Spectroscopy and QPI imaging serve as 
powerful techniques to extract the pairing symmetry. 
Meanwhile, spectroscopic imaging study of local DOS 
in the vicinity of both magnetic and non-magnetic 
impurities as well as vortices helps to distinguish 
between candidate pairing symmetries. 



Another key issue is the role of dopants. For 
most Fe-SCs, particularly the higher- 1111- and 
122-type materials, the superconductivity develops 
from the parent compounds upon chemical doping of 
electrons, holes, or even isovalent elements. In addition 
to enabling superconductivity, the dopant atoms are 
also potential sources of scattering, nanoscale phase 
separation, and electronic inhomogeneity, which may 
in turn lead to transport anisotropy, suppression, 
or vortex pinning. Therefore, STM can serve as 
an effective tool to probe these dopants on the 
nanometer scale and then understand their roles in 
superconductivity Il40ll57l . 

Finally, we anticipate that STM will play a crucial 
role in finding higher- materials. A number of 
tantalizing hints of higher-re phases, such as 53 
K Tc in single layer FeSe on a SrTiOs substrate 
[,58 J and 49 K in ~ 10% volume fraction 
of Pr-doped CaFe2As2 lf59l l60l . have appeared in 
recent months. Just as STM served to disentangle 
the phase separation and identify the true chemical 
composition of the = 32 K superconducting state 
in KFe2Se2 [3T, '44], STM should be employed to 
find the minority volume fraction superconductivity in 
rare-earth-doped CaFe2As2. Such attempts in both 
Cao.83La{).i7Fe2As2 [61 [ and Cai_i.PrvFe2As2 (where no 
clear superconducting gap is observed) have so far been 
unsuccessful, but efforts are ongoing. 

In conclusion, four years after the discovery 
of Fe-SCs, STM has contributed much to the 
understanding of these materials. Compared with 
the cuprates, Fe-SCs exhibit extraordinarily rich 
phenomenology, and can serve as a foil to unravel 
the secret of high-Tc superconductivity. There remain 
several important open questions, to which STM has the 
unique potential to provide answers. 
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